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0^ ; ABSTRACT 

Qh' We consider a class of N=l supersymmetric extensions of the Standard 

(D ! 

r-| , Model in which the soft breaking sector is CP conserving at the GUT scale. 

>► . We study the question of whether the presence of explicit CP violation in the 

^ ■ Yukawa sector of the theory induces through renormalization effects CP vio- 

' lating phases in the soft terms, which could lead to observable effects. A clear 

pattern appears in the structure of phases in the soft sector. In particular, 
the inclusion of intergenerational mixing induces large phases in the flavour 
mixing entries of the trilinear soft breaking terms, whereas the diagonal en- 
tries remain real. A mechanism is proposed for generating through chargino 
exchange a contribution to the neutron electric dipole moment which can be 
a few orders of magnitude larger that that of the Standard Model, although 
still out of reach of experimental tests. We comment on the possible relevance 
of these phases for baryogenesis at the weak scale in minimal supersymetric 
scenarios, recently considered in the literature. 



Introduction 



In this letter we reconsider the problem of CP violation in supersymmetric extensions 
of the standard Glashow-Weinberg-Salam model (SM) of the electroweak interactions. In 
particular, we want to address the question of whether the presence of a CP violating 
phase (SxAd) in the Yukawa sector of the supersymmetric theory may induce, through the 
running of the relevant parameters from the GUT scale to the Fermi scale, CP violation 
in the soft breaking sector. 

In general, it is well known that additional phases may appear in the minimal super- 
symmetric version of the SM (MSSM) due to the complexity of the soft SUSY breaking 
couplings. On the other hand, it is also well known that the presence of these extra 
phases (collectively Ssojt) would induce a 1-loop contribution to the electric dipole mo- 
ment (EDM) of the neutron, which for squarks and gluino masses of 0(100 GeV) is too 
large and requires 6soft < 0.01 [|l]. From this point of view it may be natural demand- 
ing the absence of these phases altogether by assuming that the soft breaking sector, as 
derived from the flat limit {Mpianck ^ oo) of an underlying supergravity theory, is CP 
conserving. However, the presence of an explicit CP violation in the Yukawa sector of 
the "effective" theory at the GUT scale may reintroduce some CP violation in the soft 
breaking sector at the Fermi scale, thus posing again the question of its relevance for 
present day CP odd observables. 

As a matter of fact, the presence of a small explicit CP violation Sgoft — 10^^ — 10"^ in 
the trilinear component of the soft breaking potential has been advocated by the authors 
of ref. in order to trigger the baryon-antibaryon asymmetry in the early universe in 
a finite temperature MSSM scenario. It is worth stressing that the correct sign of the 
baryon-antibaryon asymmetry depends crucially on the sign of the explicit soft breaking 
phase. It is therefore interesting to study whether the 6km induced phases previously 
mentioned have the correct size and sign to support such a scenario. 

Our main conclusion is that within the simplified approach of ref. [H, where flavour 
mixing in the effective potential is neglected, no appreciable phases are induced in the 
relevant soft breaking parameters. It turns out, however, that due to the presence of 
flavour mixing large CP violating phases may be induced in the off-diagonal components 
of the trilinear soft breaking couplings. We study the implications of such a large phases 
for the EDM of the neutron and find that they induce a contribution consistent with 
the present upper bound, although it can be as much as four orders of magnitude larger 
than the SM contribution to the elementary EDM of the quarks. The relevance of these 
off-diagonal components for the analysis of ref. deserves a separate discussion. 
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The Minimal Supersymmetric Extension of the Standard Model 



We are here interested in the class of N=l supergravity based minimal supersym- 
metric extensions of the SM, in which the spontaneous breaking of SU{2) x U{1) gauge 
symmetry is achieved via radiative corrections (RMSSM) [Q]. These models, implemented 
in a Grand Unified (GUT) scenario, exhibit the least number of free parameters and are 
consequently the most predictive ones. 

The SUSY-GUT theory is obtained from the Mpianck cxo hmit of spontaneously 
broken N=l supergravity the flat limit leaves a globally supersymmetric lagrangian, 
corresponding to the supersymmetrization of the chosen GUT model (minimally an SU (5) 
GUT) which is explicitly broken by soft terms 0]. After integration of the heavy GUT 
fields the soft terms can be casted in the following form: 

-J^^l^i^J ~ Uta{z) + Bij,hih2 - ^XaK + h.c.) (1) 

ij a 

where denote all scalar fields present in the theory. The first term is a mass term 
common to all the scalars, whereas /r is the part of the superpotential that extends the 
standard Yukawa potential 

/r = HlDlTnDl + HlEiVEEl + H^.UiVuUl 
- {H^UlVDDl + H^uiVEEl + HtDlVuUl) 

where the three r^, [x = U, D, E) are 3x3 matrices in flavour space, and the upper 
index t indicates trasposition in flavour space. Notice that in eq. (|l]) the dimensionless 
couplings r have been replaced with the massive parameters F^, and the superfields with 
their scalar counterparts. The third term of eq. (|I]) arises as the scalar counterpart of the 
term 

liH,H2 = fi{H',H', - HtH^) (3) 
present in the superpotential, while the last one is a mass term for the gauginos. 

If one assumes having a flat Kalher metric, the form of the soft supersymmetry 
breaking terms turns out to be quite simple at the scale of local supersymmetry breaking, 
that we assume to be the SUSY-GUT scale required by gauge couplings unification, 
namely Mx ~ 3 ■ 10^^ GeV . In fact at that scale we have: 

Ti = AG-T, x = D, E,U (4) 

where Aq is a massive universal coefficient (henceforth the subindex G denotes GUT scale 
quantities). In addition, each scalar in the theory gets the same mass term 

ml ■ z*Zi (5) 
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This form ensures the absence of large flavour violating effects in neutral currents at the 
low energy scale. We will also assume that the three gaugino masses are equal at this 
unification scale: 

M^ = Mg a = 1, 2, 3 (6) 

Following a widely used notation we will distinguish the superpartner fields using a 
tilde (for instance ul is the scalar partner of the left up quark ul, both belonging to the 
superfield Ul)- 

Complex Parameters in the RMSSM 

We have already mentioned that the 4 parameters Aq, Bq, Mq, fie can be a-priori 
complex parameters. In fact, by a i?— rotation, with i?— charges Qr = 1 for lepton and 
quark superfields and Qr = for the vector and the Higgs superfields, the gaugino mass 
parameter Mq can be made real; moreover, multiplying by a common phase the 2 Higgs 
superfields (with opposite hypercharge) Bc-fJ^G becomes real as well. We conclude that, in 
addition to the usual Kobayashi-Maskawa (KM) phase there are at most two tipically 
supersymmetric phases that are physically relevant, say 

arg(ylG) and arg(5G) = - arg(/iG)- (7) 

These two parameters have an important impact in the phenomenology of CP violation. 
They can induce an electric dipole moment of the quarks and leptons at the 1-loop level 
through diagrams of the type shown in fig. 1 (let us recall that in the SM the elementary 
edm of quarks arises at the third loop [0). If the masses of the particles running in the 
loop are O(IOO) GeV the EDM of the neutron is predicted to be 

rf^^^^ = 0(10-23) sin e cm (8) 

where 6soft is a typical SUSY phase (say arg(A)). This prediction has to be compared 
with that obtained in the Standard Model from the elementary EDM of the up and down 
quarks (for a review see ref. 0) : 

d^^iquarks) = 0(10-^^) e cm. (9) 

As a matter of fact, it is likely that the neutron EDM is dominated by long-distance (LD) 
effects [H, which lead to: 

dl^\LD) = O(10-=^2) e cm. (10) 
Present experiments give the upper bound 0: 

d^^^ < 12 ■ 10-26 e cm (11) 
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which, as aheady mentioned, excludes values of the supersymmetric phases larger than 
10~^ (or else requires squark masses to be at the TeV scale [0); A measure of a non-zero 
value of the EDM of the neutron in the next generation of experiments would certainly 
be a signal of new physics, among which a supersymmetric scenario with superpartners 
at the reach of the future hadron colliders. 

Running of the Soft Breaking Parameters 

If one believes that squark and gluinos are just around the corner, then a most 
conservative and natural assumption for CP violation in the soft sector is that the two 
phases in eq. (0) vanish identically at the GUT scale, due to CP conservation in the 
sector responsible for SUSY breaking. However in this case, attention must be paid to 
the KM phase, explicitly present in the model, which may induce analogous phases in the 
soft breaking sector through the renormalization of the soft breaking parameters. 

The complete set of 1-loop renormalization group equation (RGE) for the SUSY 



model here considered can be found for instance in ref. [|I^ and in explicit matrix form 



in ref. [|T2|. From inspection of the relevant RGE one easily realizes that the gaugino 
masses do not change their phase during the evolution, so that, within our hypothesis, 
they remain real; an analogous conclusion holds for the parameter /i, which depends on 
the matrices P only through the trace of the hermitian and non-negative combinations 

= ^P^ -rl X = D, E or U . (12) 

In order to discuss the evolution of the parameters B and P;^ we find convenient to 
define the matrices A^^ in the following way: 

A^ = T^- T-^ X = D, E, or U (13) 



Defining A = dA/dt, with 



t = -^-log(^) (14) 



4 



we obtain the following RGE: 

Ae = 2 (SaaMa + 3aiMi)/ 

+ 2i:i{AEaE + ^ADaD)I 

+ hUEAE + AEUE 

Au = 2 if asMs + SaaMa + f aiMi)/ 

+ 2 Tr{3Auau)I 

+ bauAu + Auau + aoAu - Auao + 2ADaD (15) 

An = 2 (f asMg + SasMa + |aiMi)/ 

+ 2 Tr(Aijaij + SAdOd)/ 

+ 5a£)A£) + + ajy^B — A/jac/ + 2Af7a;7 

B = 2 (SaaMa + aiMi) 

+ 2 Tr(A£;a£; + SAzjOd + ^Auau) 

The Aj. are a-priori generic 3x3 matrices, and the form of the initial conditions is given 
in eqs. 



Notice that one can study independently the two cases a) Mq = 0, Ac ^ and 
b) Mg 7^ 0, Ag = 0, and that the GUT value for B, namely Bq, is always additive. 

In order to analyze the evolution of the parameters we resorted to a numerical study 
of this system of RGE, coupled to the RGE for a) the gauge coupling constants ai, 
i = 1,2 or 3, b) the gaugino masses Ma and c) the "Yukawa couplings" a^, x = U,D, E. 

The initial values of the matrices at the Mz scale are calculable once the value 
of tan j3, the value of masses of the quarks and of the lepton and the Kobajashi-Maskawa 
matrix are assigned. In fact, by writing the matrices r^, in biunitary form, 

= Li^^R, (16) 

where are diagonal non-negative matrices, and making unitary redefinitions of the 
quark and lepton superfields, one finds that the parameters in eq. @) can be chosen to 
be: 

= lE 

Td = ID (17) 
Tu = K'lu 

where K is the 3x3 Kobayashi-Maskawa matrix {down- diagonal basis). 

Denoting by Vi the two vacuum expectation values of the Higgs fields, i.e. = Vi 
{i = 1,2), the matrices 'jx are related to the (diagonal) mass matrices for leptons and 
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quarks as following 

7£; = Me/Vi 

-IB = Md/vi (18) 

7(7 = Mu/V2 

This basis is very useful when performing the RGE analysis since the initial conditions 
are given in terms of masses and KM mixings and it leaves the SU{2) l symmetry explicit, 
having applied the same unitary rotation to the U and D superfields. 

To switch to the superfield basis in which the terms of the superpotential with the 
neutral Higgs superfields are flavour- diagonal {i.e. the quark mass eigenstate basis), we 
just have to redefine the Ul flavour multiplet of superfields as 

Ul K^Ul . (19) 
In this basis the mixing matrix K appears 

i) in the interactions involving the charged vector superfields 

ii) in the interactions of the charged Higgs superfields in eq. (0), 

iii) and in the analogous terms for the charged Higgs fields in eq. (|I|). 



Our numerical results will be shown in this basis. 

In order to obtain the numerical solution of the system of RGE we used the following 
set of electroweak input parameters: 



1) tan (3 = V2/V1 in the range 1 40, 



2) the MS running quark masses at mz (for a recent review see ref. [T^); in particular, 
we tested for m^p in the range 100 -r- 200 GeV. 

3) the values of the three KM angles and phase consistent with the results of the 
analysis reported in ref. |jl4|. 



At the GUT scale we spanned values for the SUSY soft breaking parameters which lead 
to a consistent SUSY mass spectrum at the electroweak scale. 

We wrote the set of the relevant RGE in matrix form and left to a program of symbolic 
manipulation the task to expand them in vector form, suited for using the variable-step 
NAG algorithm for the solution of ordinary differential equations. The results of our 
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numerical analysis revealed that, for any choice of the parameters, the amount of the 
imaginary part of B and of the diagonal elements in A^, induced at the Mz mass scale 
is completely negligible (ImjAjj}/ RejAjj} < 10~^°), whereas the off-diagonal elements of 
Au and Ad can have large phases (although the matrices Au and Ad remain, up to 1 
part per 10^^ — 10^*, hermitian). 

The following numerical example illustrates the tipical patterns we found. Consider 
for instance tan (3 = 10 and = 130 GeV. Then the values of the Yukawa matrices 
at the Mz scale read 

7C7 = diag(0. 14 10-^ 0.40 10-2,0.75) 

-fD = diag(0.29 10-^ 0.55 10-2,0.17) (20) 
7e = diag(0.29 10-^ 0.75 10-2,0.98 10-1) 



while for the KM angles and phase let us take 

i0u,023,e,,,6KM) = (0.221,0.043,0.005,0.86) 



(21) 



where we tried to maximize the CP violation effect from the KM matrix through large 
values of the mixing between the first and the third family and 6km of order 1, according 
to the present experimental constraints (see for instance ref. [|14| ). 



Finally, for Aq = 1 GeV, Mq = (notice that in this case the RGE solutions for 
scale as Aq), and Be arbitrary we obtain 



/ 



A 



u 



Ad 

Ae 
B 



\ 



(-1.3 -iO. 7)10- 
0.78, 
-1.7 10-^ 
(2.2-^1.3)10-^ 

0.99, 
(3.1-z0.1)10-3, 



(-1.0 -il. 7)10-5 \ 
-1.7 10-^ 
0.56 

;-5.1 + z3.1)10-^ \ 
(3.1 + ^0.1)10-3 
0.90 



(22) 



0.78, 
;-1.3 + ^0.7)10-^ 
V (-1.0 + a.7)10-^ 
/ 0.99, 

(2.2 + il.3)10-^ 
(-5.1 -i3.1)10-^ 
diag(0.99, 0.99, 0.98) 
-0.12 + 5g 

Increasing the top mass up to 200 GeV affects little the elements of Au, while it modifies 
by factors 2 — 3 the entries of Ad- In addition, due to the 1/ cos/5 (1/ sin/?) dependence 
of the down-quark (up-quark) Yukawa couplings, the off-diagonal elements of Au increase 
roughly as (tan/?)^, while those of Ad remain roughly constant. The two features above 
are in fact related and may be qualitatively understood by analizing the interplay between 
up an down Yukawa couplings in the RGE for Au and Ad, (eqs. (|T5|)). 

As mentioned before, the matrices Au and Ad turn out to be quasi- hermitian (non- 
hermiticity appears at the 1 part per 10-'^ — 10-^ level in the off diagonal elements and we 
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have neglected it). Notice also that the matrix Ae remains in any case diagonal, because 
of the absence of any source of lepton flavour violation. 



A consequence of the "hermiticity" of the matrices is the reality of B (up to a 
part in 10^^°). In fact, the evolution of B depends on the trace of the product of two 
hermitian matrices {ay ■ Ax), which is real. 

The fact that the diagonal entries of A^ are to an extremely good approximation real 
bears important implications for the phenomenology of SUSY induced CP violation, which 
we will shortly discuss. This feature can be understood as follows. Suppose decomposing 
the matrix A^ in its hermitian and antihermitian parts A^ = A'^ + A'^ respectively. It 
is then evident that the boundary conditions for A^ at the GUT scale, eq. (H), imply 
A"q = 0, so that the antihermitian part can be generated only through [A'^, but, due 
to the diagonality of A'^ q also the commutator is zero at the GUT scale. Moreover, to get 
a non-zero value of {A")ii from this commutator one readily verifies that a certain amount 
of non- diagonality in {A'^)ij, i 7^ j, is also needed, which is only induced radiatively. 

The Neutron Electric Dipole Moment 

We shall now consider some consequences of the previous results. Since both fi and 
B remain real, also the mass parameter mg = —Bfi responsible for the mixing of the 
two scalar higgses remains real at the end of the running. Analogously the gaugino mass 
matrices, i.e. the chargino and the neutralino mass matrices, are real; this fact implies 
that only real orthogonal matrices are needed to reach the basis of mass eigenstates. Also 
the sleptons mass matrices are real (and diagonal), so that the CP violating effects in 
the leptonic sector, such as the electron or muon EDMs are zero at 1-loop level in the 
RMSSM as in the SM. 

Let us now consider the role the complex off-diagonal parameters in the matrices A^ 
for the squark mass matrices, and begin our discussion with the up-squark matrix. The 
6x6 scalar quark mass matrix can be written in terms of 3 x 3 submatrices as follows: 

A/f2 _ I ^i,LL ^i,LR 1 /r,o\ 

The two LL and RR blocks are determined evolving the scalar mass parameters in eq. 
(|]) with the appropriate RGE; the parameters A^ enter only the LR (and RV) blocks 
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The basis of squark mass eigenstates is reached using the unitary rotation Su, defined by: 

S^MiSi = diag(mij (25) 

It is useful to spht Su into two 6x3 submatrices 

Su = {Su^L, Su^ji) (26) 

which relate the the scalar partners of the left and right-handed quarks to the scalar mass 
eigenstates. 

Notice that in the case of the down squarks, since M'-^^^ is proportional to M^, 
the unitary matrix is block diagonal up to 0{MD/msusY)', as we will see this fact 
suppresses the up-quark contribution to the electric dipole moment of the neutron. 

Let us finally come to the neutron EDM. Using a non-relativistic quark model, the 
elemntary EDM of the quarks and d^. are related to the dipole of the neutron according 
to: 

4 1 

dn = -dd - -du- (27) 

Let us consider first the d quark. The dipole can be computed from the imaginary part 
of the amplitude A^"^ for the chirality- flipping dn ^ di + J process according to 

dd = -Im(41?) (28) 

In the SUSY model there are no 1-loop contributions to the neutron EDM induced 
by or exchange. The only possibility left is to examine typically SUSY diagrams 
with squarks running in the loop (together with gluinos, charginos or neutralinos) , and 
resort to that component of the amplitude in which the needed helicity flip is realized 
in the loop. From a qualitative mass insertion analysis of the relevant diagrams it is 
easy to convince oneself that a non vanishing (and possibly complex) LR mass insertion 
in the squark line is needed. In fact, when gluinos (or neutralinos) are considered, the 
down-quark EDM turns out to be proportional at the leading order to Im{(74£))ii}, which 
is zero in the scenario here considered. One might circumvent this problem by invoking 
flavour changing (FC) effects at the gluino-quark-squark vertex so to involve complex 
off-diagonal entries in the soft insertion. However these FC effects are of radiative origin, 
and for squark and gluino masses heavyer than 100 GeV turn out to give quite a small 



contribution to the dipole amplitude (see ref. [|T2 



We are therefore lead to consider chargino exchange, with up-squarks running in the 
loop. Analogously to the gluino case if one considers the "diagonal" exchange of Uu,c,t 
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no effect arises due to tlie reality of tlie diagonal entries of tlie trilinear soft breaking 
parameters. On the other hand, there exist now also "non-diagonal" components of 
the amplitude at the leading order (see fig. 1) which are proportional to off-diagonal 
elements of the Au matrix and to complex combinations of the relevant KM mixings. In 
particular we may expect that the dominant contribution arises through the exchange of 
the "higgsino" component of the chargino field, via the quark-squark "flavour- chain" 

d-R {ul, cl) ^iji^ di (29) 

since we can take advantage of the presence of the large top Yukawa coupling. For large 
tan/?, then, also the presence of the down Yukawa coupling becomes important as an 
enhancement factor. 

More concretely, using the interaction lagrangian (see for instance the appendix of ref. 

H) _ 

+ PniUl^S^^LK^d)] d + h.c. 
where U and V are the 2x2 orthogonal matrices responsible for the diagonalization of 
the chargino {'XT) mass matrix (see for instance ref. ||l5l), we find that the contribute to 
the d-quark EDM is can be written as 



The function T[x) is given by 



^hu [sin:F K^d 



e cm (31) 
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J^(x) = _J__(5 - 12a; + 7x^ + 2x{2 - 3x) ln(x)) (32) 
6(1 — x)'^ 



A formula analogous to eq. (|3TD holds for d^. However, the appearance of the down 
quark mass matrix in AfJ^^ (compare with eq. (|2^) ) suppresses the EDM of the up-quark 
compared to that of the d-quark. 

To numerically estimate this effect we considered the lightest chargino to be close to 
the present experimental limit, say m^^ ~ 50 GeV, lightest squarks masses of 0(100 GeV) 
and maximize the amplitude in the remnant parameters. We therefore find that the SUSY 
contribution to the neutron EDM in this class of models generated by the elementary EDM 
of the quarks can be as large as 

^susY ^ 0(10-30) ('-^\ e cm (33) 
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where the hnear dependence on tan (3 for large (3 comes from the presence of the down 
quark Yukawa couphng (the dependence on the top mass in the studied range is weaker). 

This result has to be compared with the analogous SM contribution of eq. @. 
In spite of the absence of phases in the diagonal entries of the trilinear soft breaking 
terms, the smaller off-diagonal terms may induce an elementary quark EDM four orders 
of magnitude larger than that of the SM, and give a neutron EDM still two orders of 
magnitude larger than that induced by long distance effects, eq. (^). It is obvious, 
however, that this mechanism cannot explain a measurement of the neutron EDM close 
to the present experimental bound. 

For what concerns the possible implications of these results for the analysis of ref. 
1^ , on the generation of baryon-antibaryon asymmetry at the weak scale in a finite tem- 
perature MSSM scenario, no phases of 0(10^^ — 10"^) are induced in the diagonal (top) 
entries of the trilinear terms as suggested by the authors. However, the aforementioned 
analysis just neglects intergenerational mixing in the calculation of the effective potential. 
The possibility that the large phases induced in the off-diagonal entries of may actu- 
ally become relevant and trigger the transition to a baryon dominated universe remains 
open. Answering this question requires a detailed analysis which is beyond the scope of 
the present letter. 

Figure Captions 

Figure 1. The leading SUSY contribution to the elementary EDM of the quarks in the 
class of models described in the paper is shown in the interaction eigenstate basis. The 
photon is attached in all possible ways. 
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